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ABSTRACT 

In  this  paper  we  examine  two  couplings  that  produce 
time-shifted  synchronization  in  a  pair  of  chaotic 
oscillators.  The  couplings  do  not  utilize  an  explicit  time- 
delay  term.  We  characterize  the  observed  synchronization 
by  determining  the  dependence  of  the  time  shift  and  cross 
correlation  between  the  drive  and  response  oscillators  on 
a  tunable  parameter.  Our  observations  agree  well  with 
estimates  of  the  time  shift  and  cross  correlation  using  a 
transfer  function. 


1.  INTRODUCTION 

Coupled  chaotic  oscillators  often  display  some  form 
of  synchronous  behavior  [1,  2].  Recently,  varieties  of 
synchronization  have  been  identified  where  the  oscillator 
waveforms  are  nearly  but  not  exactly  identical  [3-9].  In 
these  cases,  either  the  coupling  or  a  parameter  mismatch 
produces  a  time-shifted  synchronization  state  with  some 
finite  amplitude  error  (e.  g.,  see  Fig.  1).  It  remains  an 
interesting  open  question  whether  natural  systems  utilize 
such  mechanisms.  On  the  other  hand,  in  engineered 
systems  applications  have  been  identified  where  an  easily 
varied  time  delay  is  desired. 

Recently,  chaotic  oscillators  have  been  suggested  as 
efficient  sources  for  generating  wide-bandwidth 
waveforms  for  various  applications.  For  example,  the 
broadband  and  non-repeating  nature  of  chaos  provides  an 


Figure  1.  Example  of  time-shifted  synchronization  in  a 
pair  of  coupled  Rossler  oscillators.  The  response 
waveform  ^/(t)  lags  the  drive  waveform  X](t)  and  is 
slightly  attenuated. 


ideal  combination  of  high  range  resolution  and  no  range 
ambiguity  for  radar.  Local  coupling  can  synchronize  an 
extended  array  of  chaotic  oscillators,  thereby  providing  a 
coherent  state  suitable  for  power  combining  in  a  wide- 
bandwidth  phased  array  antenna.  Moreover,  time-shifted 
synchronization  provides  a  mechanism  to  steer  the 
radiated  beam.  A  tunable  coupling  strategy  may  be  used 
that  continuously  varies  the  time  shift,  thereby 
electronically  steering  the  beam  to  a  desired  direction  [3, 
5,  6,  9].  Effectively,  these  systems  trade  some  decrease  in 
cross-correlation  between  oscillators  for  a  practical  means 
of  obtaining  a  tunable  time  shift.  In  practice,  this  tradeoff 
appears  as  reduced  beam  quality  for  larger  steering  angles 
and  limits  the  maximum  steering  angle  for  these  arrays. 

In  this  paper,  we  examine  two  unidirectional 
coupling  schemes  that  produce  non-trivial  time  shifts  and 
cross-correlations,  and  we  compare  these  observations  to 
predictions  based  on  a  new  method  for  estimating 
characteristics  of  the  synchronized  state  [10].  This  method 
uses  a  transfer  function  derived  from  the  coupling  model 
to  provide  estimates  of  the  time  shift  and  amplification  or 
attenuation  between  the  drive  and  response  waveforms. 
Each  coupling  scheme  involves  a  parameter  to  tune  the 
time  shift,  and  typically  a  large  range  of  time  shifts  can  be 
obtained.  Overall,  we  find  that  observed  time  shifts  and 
cross-correlations  agree  well  with  estimates  based  on  the 
transfer  function  and  that  the  agreement  increases  with 
coupling  strength. 

2.  TIME-SHIFTED  SYNCHRONIZATION 

We  consider  a  pair  of  identical  chaotic  oscillators 
coupled  unidirectionally  of  the  form 

x  =  f(x),  y  =  f{y)  +  g(x,y )  (1) 

where  x(t)  and  y(t)  are  the  state  vectors  of  the  drive  and 
response  oscillator,  respectively, /is  a  chaotic  flow,  and  g 
is  a  linear  coupling  function.  Since  /  is  chaotic,  it  is 
necessarily  nonlinear  and  the  state  vectors  x  and  y  each 
contain  at  least  three  scalar  states. 

Here,  we  consider  g  of  the  form 

g(x,y)  =  k{Pxl-Qyl}el,  (2) 

where  x,  and  y\  are  the  first  components  of  the  state 
vectors  x  and  y,  respectively,  k  is  the  coupling  strength, 
et  =  (l,  0,  0)r ,  and  P  and  Q  are  linear  scalar  operators. 

At  least  five  coupling  schemes  of  this  form  have 
previously  been  identified  that  produce  approximate  time- 
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shifted  synchronization  [6-8,  10  ,11].  In  these  cases, 
identical  synchronization  is  not  an  exact  solution  of  Eqs. 
(1)  and  (2).  Instead,  we  consider  the  drive  and  response 
synchronized  when 

k{Pxl  -Qyi}=s(t)  (3) 

and  s(t)  is  small  compared  to  the  first  component  of  the 


vector  field  / .  In  other  words,  the  coupling  term  only 
provides  a  small  perturbation  to  the  natural  dynamics  of 
the  response  oscillator. 

Apart  from  a  time  shift,  the  waveforms  x t  (t)  and 


y1  (t)  are  typically  similar  in  shape,  with  the  primary 

distortion  being  just  a  change  in  amplitude.  Thus,  we  can 
write 


Ti(0~  Ax,(t-r) 


(4) 


where  r  is  the  time  shift  and  A  is  the  approximate  gain 
(or  attenuation)  of  the  response  waveform  relative  to  the 
drive.  For  a  given  lag  St ,  we  define  the  cross  correlation 
function 


Jvj  (t)yl  (t  +  St)dt 

J  v,  (t)vj  [t)dt 


(5) 


where  we  note  the  normalization  does  not  treat  the  drive 
and  response  waveforms  symmetrically.  Assuming  Eq. 
(4),  we  recognize  8t  =  x  maximizes  the  cross  correlation 
and 


C(t)«A 


(6) 


so  that  the  cross  correlation  provides  an  estimate  of  the 
response  gain. 

It  has  been  shown  that  a  reasonable  estimate  of  the 


time  shift  r  and  cross  correlation  C(r)  can  be  derived 


by  neglecting  the  small  perturbation  s(t )  [10].  Thus, 


setting  s(t)  =  0  in  Eq.  (3)  and  taking  a  Fourier  transform, 
we  can  write 


(a)  (b) 


X 


Figure  2.  (a)  Correlation  and  (b)  time  shift  between 
coupled  Rossler  oscillators  observed  in  numerical 
solutions  (circles)  and  estimated  using  the  transfer 
function  Eq.  (12)  (line). 


j)1(®)  =  r(®)i1(®)  (7) 

where  (®)  and  jq  (a>)  are  the  Fourier  transforms  of 
X[(t)  and  y,(t),  respectively.  We  note  that  the  transfer 
function  r(©)  depends  only  on  the  linear  operators 
P  and  Q  in  this  approximation,  and  that  for  many  simple 
operators  T  (®)  can  be  calculated  analytically.  The  time 

shift  between  the  drive  and  response  is  then  predicted  to 
be 

rArgFK))  (8) 

®o 

while  the  predicted  cross  correlation  is 

C(r)  =  |2>0)|  (9) 

where  ®0  is  the  dominant  angular  frequency  of  the 
specific  chaotic  oscillator  under  consideration,  and 
arg(T(®0))  represents  the  phase  angle  of  the  complex 

transfer  function  evaluated  at  this  frequency. 

Below  we  introduce  two  examples  of  couplings  that 
produce  approximate  time  shifted  synchronization.  Each 
coupling  includes  a  tunable  parameter  that  varies  the  time 
shift  and  cross  correlation  of  the  synchronized  state.  We 
determine  the  functional  dependence  of  these  quantities 
on  the  tunable  parameter  and  compare  numerical  results 
with  predictions  made  using  a  transfer  function. 


3.  FIRST  EXAMPLE 


We  first  consider  a  pair  of  Rossler  oscillators  whose 
flow  is  determined  by  the  vector  field 


/(*)  = 


x,  +  0.15x2 
0.2  +  (x,  -10 )x3y 


and  who  are  coupled  by  the  function 

g (x,  y)  =  k  {*i  (?)  -  y,  (t)  +  Ayt  (t )}  e, , 


(10) 


(11) 


where  k  =  0.3  and  A  is  a  tunable  parameter.  Typical 
time  series  of  xx  and  )\  with  A  =  -0.5  ,  shown  in  Fig.  1, 
exhibit  a  state  of  time-shifted  synchronization.  The 
primary  form  of  distortion  of  the  response  trajectory 
relative  to  the  drive  is  a  small  degree  of  attenuation. 

Following  the  procedure  described  in  the  preceeding 
section,  we  obtain  the  transfer  function 


T(co) 


Ti  H 
A  H 


1 

1-icoA 


(12) 


The  spectral  content  of  x,  (/)  is  dominated  by  the  single 
frequency  (on  »  1 .04 .  The  magnitude  of  the  transfer 
function  at  a>0 ,  given  by 


2 


|7>o)  =■ 


1 1  +  co0  A 


suggests  that  the  amplitude  of  the  response  oscillation  will 
be  attenuated  relative  to  the  drive  when  |/l|  ^  0  .  The 
phase  delay  of  the  transfer  function,  divided  by  ®0 
predicts  the  time  shift  of  the  synchronized  state.  In  this 


arg(^(®0))  1  t  ^ 

- 1 - -  =  — tan  {(o0A), 


Notably,  at  large  values  of  \A\  the  predicted  time  shift 
asymptotes  to  ±/r/(2®0) . 

We  compare  these  theoretical  estimates  to 
observations  of  the  correlation  function,  Eq.  (5).  The 
value  of  St  at  which  C[St)  reaches  its  first  maximum  is 
identified  as  the  observed  time  shift  r  .  The  magnitude  of 
C(<?t)  at  its  maximum  is  a  measure  of  the  gain  or 

attenuation  of  the  response  waveform  relative  to  the  drive. 
The  time  series  data  is  generated  by  numerical  solution  of 
Eqs.  (1),  (10),  and  (11). 

Fig.  2(a)  shows  the  magnitude  of  the  transfer 
function,  T(®0  )|  (solid  line),  agrees  well  with  the 

observed  correlation,  C  ( r )  (circles),  over  a  wide  range 
of  values  of  the  parameter  A  .  Interestingly,  the  agreement 
is  good  even  for  large  values  of  A  where  C(r)<0.4. 
This  is  surprising  since,  under  these  conditions,  the 
approximation  that  ^(t)  is  negligibly  small  is 
presumably  rather  poor. 

The  predicted  values  of  the  time  shift  also  agree  quite 
well  with  the  observed  time  shift  r ,  as  shown  in 
Fig.  2(b).  Most  important,  the  observed  asymptotic  levels 
of  the  time  shift  agree  with  the  values  predicted  by  Eq. 
(14)  to  within  a  few  percent.  Again,  it  is  surprising  that 
the  agreement  is  so  good  at  large  values  of  \A\ . 

Thus  far  we  have  considered  only  a  single  value  of 
the  coupling  strength,  k,  and  have  not  addressed  the  role 
of  coupling  strength  in  determining  the  properties  of  the 
synchronized  state.  Obviously,  the  coupling  strength  must 
be  sufficient  to  guarantee  stability  or  the  properties  of  the 
synchronized  state  will  not  be  observable  [1,2].  However, 
there  typically  exists  a  range  of  coupling  strengths  over 
which  synchronization  is  stable,  and  the  question  remains 
in  this  case  whether  the  coupling  strength  further  affects 
the  time  shift  and  cross  correlation.  One  might  guess  that 
stronger  coupling  would  tend  to  minimize  the  difference 
term  Pxl  -  Qyx  in  Eq.  (1).  Thus,  the  properties  of  the 
response  state  y,  would  more  closely  agree  with  the 
transfer  function.  Here  we  examine  a  particular  example 
that  is  consistent  with  this  hypothesis.  (We  note  that  the 
effects  of  increased  coupling  on  the  other  states  of  the 
response  may  be  quite  different.  In  fact,  the  resemblance 


Coupling  Strength 


Figure  3.  Root  mean  square  error  in  the  transfer 
function  predictions  of  the  (a)  correlation  and  (b)  time 
shift  as  a  percentage  of  the  observed  value. 

of  y2  and  to  their  drive  counterparts  may  actually 
decrease  with  increasing  coupling  strength.) 

We  focus  on  the  range  of  coupling  strengths  from 
k  =  0.1  to  k  =  2.0.  The  lower  bound  on  this  range  is 
determined  by  the  requirement  that  synchronization  must 
be  stable.  We  find  that  the  transition  between 
synchronized  and  unsynchronized  dynamics  occurs 
somewhere  around  k  =  0.1  .The  upper  bound  is  set  by  the 
fact  that  when  k  -  A'1  the  derivative  term  in  the  coupling 
exactly  balances  the  derivative  term  on  the  left  hand  side 
of  the  equation.  This  cancellation  changes  the  character  of 
the  equations  and  causes  problems  for  numerical 
simulation.  Thus  we  avoid  it  by  constraining  k  to  fall 
below  A~'  for  the  particular  value  of  A  we  wish  to 
consider.  We  fix  the  value  A  =  -0.5  under  the  assumption 
that  typical  behavior  occurs  at  this  value. 

Following  the  same  procedures  as  described  above, 
we  can  numerically  integrate  this  system  and  observe  the 
time  shift  and  cross  correlation  as  the  coupling  strength 
k  is  varied.  We  quantify  the  agreement  between  the 
observed  values  and  the  predictions  of  the  transfer 
function  by  calculating  the  root  mean  square  deviation 
between  the  observed  and  predicted  values  as  a  percent  of 


3 


the  observed  value.  Fig.  3  shows  the  prediction  error  in 
the  (a)  correlation  and  (b)  time  shift  as  functions  of  the 
coupling  strength.  As  expected  the  prediction  error 
decreases  with  increasing  coupling  strength  over  the 
entire  range.  Thus  the  transfer  function  provides  a  better 
prediction  as  the  coupling  strength  is  increased. 


4.  SECOND  EXAMPLE 


We  next  examine  a  variation  of  the  previous  coupling 
scheme  that  includes  a  second  derivative.  Specifically,  the 
coupling  scheme  is 

g(x,y)  =  k(xl-yl+Ay,+ayl)el,  (15) 


where  A  is  a  tunable  parameter  and  a  is  a  fixed 
parameter.  The  class  of  oscillators  introduced  by  Sprott  is 
one  for  which  this  coupling  is  easily  implemented  [12]. 
Here  we  use  one  such  oscillator  whose  flow  is  defined  by 


By 


f(x)  = 


x3 


(16) 


0.6X3  ~X2  +  |*i|-2 j 

numerical  solution  we  determine  the  dominant 


frequency  in  the  spectral  content  of  x,  (/)  to  be 


0.958  ±  0.005  .  The  transfer  function  for  this  coupling  is 
T(co)  =  \-aco2  +icoA.  (17) 

The  predicted  cross-correlation  is 

C(r)  =  ^{\-aco2)2 +co2  A2 .  (18) 


Notably,  the  correlation  at  A  =  0  is  predicted  to  be 
1  -  aa>02 .  All  previous  couplings  that  have  been  analyzed 
using  the  transfer  function  produced  identical 
synchronization  at  A  =  0  [6-8,  10  ,11].  The  predicted 
time  shift  is 


1  _! 
r  =  — tan 
CO n 


co0A 
1  -  acon 


(19) 


(a)  (b) 


Figure  5.  Example  of  time-shifted  synchronization  in  a 
pair  of  coupled  Sprott  oscillators.  The  response 
waveform  jy(t)  anticipates  the  drive  waveform  Xy(f). 

We  numerically  integrate  the  system  to  generate  time 
series  of  x,  and  y,  at  each  of  a  range  of  values  for  A 
and  with  k  =  20  and  a  =  0.5  .  From  the  numerical 
simulation,  we  estimate  the  time  shift  and  attenuation 
from  the  cross-correlation  function,  Eq.  (5).  Fig.  4  shows 
the  observed  dependence  (circles)  of  (a)  C(y)  and  (b)  r 
on  the  parameter  A  .  The  observed  time  shift  dependence 
on  A  follows  closely  the  prediction  of  Eq.  (19).  The 
observed  cross-correlation  at  A  =  0  is  0.56  ±  0.005  which 
agrees  well  with  the  prediction  of 
l-aa>02  =0.541  ±0.005 .  The  correlation  then  grows 
quadratically  as  \A\  increases,  in  agreement  with  Eq.  (18) 
and  reaches  unity  around  |d|  =  0.9  .  Fig.  5  shows  a  typical 
time  series  of  x,  and  y,  with  A  =  0.9  which  verifies  the 
absence  of  significant  distortion  despite  the  time  shift. 


5.  CONCLUCIONS 

We  have  shown  two  examples  of  couplings  that 
produce  time  shifted  synchronization  and  verified  the 
accuracy  of  the  transfer  function  in  predicting  the 
observed  properties  of  the  synchronized  state.  These 
results  represent  further  evidence  of  both  the  wide  range 
of  coupling  that  can  produce  time-shifted  synchronization 
and  the  usefulness  of  the  transfer  function  in  analyzing 
this  phenomenon.  In  the  first  example,  we  observed  that 
the  accuracy  improves  with  increasing  coupling  strength. 
This  result  suggests  that  strong  coupling  would  be 
preferable  in  applications  where  that  the  transfer  function 
is  utilized  as  a  design  tool  [10]. 
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